close all; clear; clc;

I=imread('lina.bmp');
F=fft2(I);     %频谱图
Fc=fftshift(F);  %中心平移频谱
Fc1=fftshift(F);
[r,c]=size(Fc);
Fc1(floor(r/2-59:r/2+60),floor(c/2-59:c/2+60))=0;
%在中心平移后频谱图中截取中心120*120区域
iFc1=ifft2(Fc1); %反变换效果图
subplot(2,2,1);imshow(I);title('原图像'); 
subplot(2,2,2);imshow(Fc,[]);title('中心平移频谱图');   %[]增加对比度
subplot(2,2,3);imshow(Fc1,[]);title('去低频频谱图');
subplot(2,2,4);imshow(abs(iFc1),[]);title('反变换效果图');